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ABSTRACT 


Nowadays, wind turbine energy has an increased importance in electrical 
power applications since when it is considered as an essential inexhaustible 
and broadly available energy resource. An aerogenerator is a device that 
transforms a part of the kinetic energy of the wind into available mechanical 
energy on a transmission shaft, and then into electrical energy through a 
generator, which is in our case a dual power asynchronous machine. 
An important characteristic of a wind turbine is that the avail, able maximum 
power is provided only in a single given operating point, called Maximum 
Power Point. Many classical methods and controllers have been widely 
developed and implemented to track the maximum power point. Among 
drawbacks of a classical PI controller is that its parameters are not constant, 
these conventional control laws may be are insufficient because they are not 
robust, especially when the accuracy requirements and other dynamic 
characteristics of the system are strict. The new idea in this paper is to 
introduce the Genetic Algorithms theory into the controlstrategy that used 


inthe conversion chain of the wind turbine, in order to improve stability. 
Simulation results approve that the application of Genetic Algorithms to the 
PI regulator, minimize or eliminate the drawbacks of the classical PI 
regulator, and greatly promote the efficiency and stability of the wind 
turbine systems. 
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1, INTRODUCTION 

For several years,traditional power generation plants cause various pollutions and ecological 
imbalance by using fossil and natural sources such as uranium, hydrocarbons, and water.For example, 
thermal power plants (coal, oil) are responsible for atmospheric emissions. Hydroelectric power stations 
sometimes require the construction of water reservoirs that radically alter the equilibrium of an entire 
area [1]. The degradation of air quality is considered as a major problem which has led authorities around the 
world to take measures to reduce pollutant emissions. As a result, the power generation industry is 
undertaking to adapt to these new constraints by promoting wind energy. Currently, several wind farms are 
installed in the world and are the only one that is both inexpensive and capable of mass production, 
Wind energy has a higher energy potential and becomes the first source of renewable energy after 
hydraulics [2]. Most variable speed wind turbine systems are equipped by the dual-powered asynchronous 
machine, this type of machine allows a wide range of speed variation (+30% ) of the synchronism 
speed [3, 4]. This is very useful to extract the maximum power during its operation. 
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Wind turbines have a specific electrical characteristic of a non-linear nature and have a particular 
point called "Maximum Power Point". The purpose of any control in this domain is the permanent search for 
that’s maximum, this technique is also widely used in photovoltaic systems [5, 6], where the MPPT's search 
in the current-voltage plan ultimately allows the cost of the kWh produced to be reduced. This power 
optimization occurs in zone 2 of the operating curve of a variable speed wind turbine. In recent decades, 
several advanced control methods have been developed to control the production processes of electrical 
energy. However, the performance of the control system has not always been assured due to the complexity 
of the control algorithms and the non-linearity characteristic of the systems to be controlled. 

The Conventional PI control techniques cover a wide range of industrial applications; those 
techniques are linear and have the advantage of simple implementation and the non-complexity synthesis of 
regulators’ gains. Inthe course of time, these applications have become less efficient, especially if the 
processes to be controlled have acomplex and non-linear structure. Finding an alternative to these 
conventional techniques is a persistent concern for many researchers because it involves a compromise 
between the robustness on one side and cost on the other. According to reference [7], classical PI controllers 
have some problems related to non-linear systems with non-constant parameters; the conventional control 
laws may be insufficient because they are not robust, especially when the requirements of accuracy and other 
dynamic characteristics of the systemto be controlled are strict. For this, we should use insensitive control 
laws to the disturbances and nonlinear cases. The main problem in this article is the calculation of the 
parameters of the classical PI regulator. 

There is a call in the literature that encourages the improvement of crossover operators for more 
efficient optimization schemes of evolutionary algorithms. The influence of crossing operators is vital for the 
whole process of the genetic algorithm in the search for an optimal search space [8, 9]. The novel idea of this 
article is to introduce the theory of genetic algorithms in the control strategy used in the wind turbine 
conversion chain to improve stabilityby itsapplying on the conventional PlI.wherein our casethe genetic 
algorithm is used to search for an optimum PI regulatory gains.The simulation results confirm that the 
application of genetic algorithms on the PI regulatorseliminatestheirdisadvantages and greatlyfavours the 
efficiency and stability of wind turbines. 


2. WIND ENERGY CONVERSION CHAIN MODELING 

The overall wind turbine diagram connected to the grid is illustrated by Figure 1. The wind turbine 
captures the kinetic energy of the wind and converts it into a torque that turns the rotor blades. There are 
three factors determine the relationship between wind energy and the recovered mechanical energy by the 
rotor: the air density, the area swept by the rotor and the wind speed. The Air density and wind speed are 
climatological parameters that depend on the site [10]. 
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Figure 1. The control strategy principle 


2.1. Turbine Modeling 
The available wind power crossing a surface S is defined by [11-12]: 


1 3 
ev 
(1) 


Indonesian J Elec Eng & Comp Sci, Vol. 18, No. 1, April 2020 : 167 - 178 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 Oj 169 


The available mechanical power depending on the generator speed for different wind speeds 
[13-17]: 


1 3 
Fin = 5 Cp Ap syj 


(2) 


For the turbine that used in this study, the power coefficient 1s approximated by the following 
formula [18]: 


116 


—21 
C,(A, B) = 0.5 (= — 0.48 —5) exp(=—) + 0.0068 (3) 
The speed ratio defined by: 


pe Ovurbine R 
v (4) 


The dynamic fundamental equation makes it possible to determine the mechanical speed evolution 
of the asynchronous Generator. The simplified model of this equation is given by: 


dQ 
——HHEe — Cg — Cem — Cvis 
dt (5) 


J 


3. OPTIMIZATION by GENETIC ALGORITHMS 

The Optimization is the search for the best solution to a problem in the sense of one or more 
selected criteria, respecting the characteristics of the system and the imposed constraints on it. The optimal 
solution is sought from a population of solutions by using random processes. The search for the best solution 
1s carried out based on the creation of a new solutions generation, by the application of three operators 
successively; to the current population: the selection, the crossing, and the mutation. These operations are 
repeated until a stop criterion 1s reached [19]. A simple Genetic Algorithm principle is shown in Figure 2. 












Generation of | 
the wutial 
population 





Figure 2. The general principle of genetic algorithm 


The different stages of the simple Genetic Algorithm: 
a) The initial population Creation. 
b) The individual's Evaluation of the population. 
c) The best individuals Selection. 
d) Reproduction (crossing and mutation). 
e) Anew generation Formation. 
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4. GA APPLICATION to the PI REGULATOR OPTIMIZATION, for a WIND TURBINE 
CONTROL 

In this part, we detail the implementation aspects of the PI parameters optimization by using GA. 
The aim is to find all the optimized parameters P = [KPopt, KIopt]. So that the system speed response of the 
closed-loop wind turbine is stable and as robust as possible for any variation of J, belonging to the gap. 
By using GA, it is possible to search for the optimal solution in a three-dimensional space, to the detriment of 
the quasi-one-dimensional search, which has been implemented from the iterative approach. Each one of the 
three-dimensional search space corresponds to the individual x variable of that we seek to optimize. 
The principle of the PI corrector optimization by using genetic algorithm as shown in Figure 3. 


[Xin Xin] = [Kpp Ky] (6) 





Figure 3.The principle of the PI corrector optimization by using genetic algorithm 


In the same way, that we limited the uncertainty of the mechanical load inertia Jt according to the 
defined interval in the work, a region of space D limits the possible parameters values of the regulator PI, 
defined by the best coefficients for the Jmin terminals on the one hand and Jmax on the other hand. 
From some knowledge of the system to be regulated and the design of a PI type regulator, we can limit the 
search region of the optimal solution while ensuring the approach optimization, so that the GA convergence 
time is greatly reduced. 


D=[x.. ]*[x.. J=[k .k }*[k. . ok, ] 
lin lin pmin pmax imm imax (7) 


The approach optimization of the PI regulator by GA 1s defined in the following steps: 


4.1. Initial Population Definition 
The initial population is randomly selected in a solution search space. The variation interval of each 
variable of D is defined as follows: 


4.1.1 We Adopt The End Of The Variation Range Of Jl As A Reference For The System: 
t ‘ min (1) 


J. 
t imax (2) 


4.1.2 For each reference of inertia 
we calculate from the control system simulations of a complete closed-loop wind _ turbine, 

the maximum values of each coefficient KP and KI by respecting the criterion of robustness and stability. 
As aresult, we have two sets of controller parameters: 

a) Together (1): 

) g C)ik mena 
b) Together (2): 

) g (2):K ks 


4.1.3 The search space D is defined by: 


Indonesian J Elec Eng & Comp Sci, Vol. 18, No. 1, April 2020 : 167 - 178 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 Oj 17] 


The optimal values of the PI regulator are included in the search space D: 


opt — K popt = lomax K pK po) 


‘opt — K opt - lomax Ky K7)| (9) 


4.2. The Adaptation Function 


@ 
The integral inverse of the quadratic error between the reference speed of the mechanical load "& 
and the real speed @ 1s used as an adaptation function to maximize over an interval T: 


1 
Jad Pop 


a 2 

PO @) dt 
0 (10) 

1 
Pr ay a 
Imax 12 = 
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Where L-1 represents the inverse transformation of Laplace, and * represents a temporal convolution. 
The uncertain transfer function between the motor torque and the load speed Hl-c (p, n) is 
given by: 


1 











a (p,n) = 
: is he a Ry 
el Nee ee! Us ee op +H] 
fi KG +f) BONG a ad K(f +f) 
m “1 1=m “1 l1=m “1 


(12) 


4.3. The Selection Operator 
We use the standard geometric selection method which is a rank selection method based on a 
standardized geometric distribution. This method limits the probability of the best individual selection at the 


th 
first generations beginning. With this method, the selection probability of the i 


nN 
The individual in the population of size ? is given by: 


rool 
oer 





Ps 


"p 
1-d-q,) (13) 


The rank of an individual is the index of the population in which it has been noticed as being 
dominated by no other individual. From the individuals set of a population, you take out those who are not 
dominated by any other individual, and you assign them to rank 1.On the individuals who did not stand out, 
you repeat the operation, looking for those who are not dominated by any remaining individuals and you 
assign them to rank2, and so on until the exhaustion of population. The best-judged individuals are those of 
the lowest rank. 


4.4. The Crossing Operator 
As part of our study, we use the arithmetic crossing. The crosses are considered with two individuals 


(Pap Paz) 


from the same population , generally called fathers. As result of this crossing, we have two new 


individuals (€,,€,) , which are called sons. This type of crossing is defined as follows: 


(é = Pal, + Par) 
Cy = Py A-1)+ Park (14) 
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Where © is a random weighting variable between 0 and 1.The number of times that the arithmetic crossing 


must be applied to each generation is defined by the probability of crossing re 


4.5. The Mutation Operator 
We use a real mutation with uniform type. Each individual of any population is changed according 


to a certain probability Pm by a random value, which is drawn following a uniform distribution over the 


interval of possible solutions D or the specific interval of mutation ™. 


4.6. The Iteration 
The steps of the Genetic Algorithm (4.1) - (4.4) are executed until the moment while the adaptation 
function value of individuals no longer evolving, or when the maximum number of generations 


n 
§ is reached. 


5. WIND TURBINE SIMULATION 

All simulations were performed under MATLAB / SIMULINK software environment on a wind 
system. Two MPPT techniques were implemented in the simulation as controller; to study and compare the 
dynamic response of the wind system in terms of the variable wind speed evolution. These two control 
structures were simulated by considering an average wind profile around (12.5 m/s) as shown in Figure 4. 
We display the obtained results for the different control strategies used. The power coefficient performance of 
the turbine and the relative speed (A), the generator rotation speed as well as its reference and the turbine 
aerodynamic power. The proposed (PI-genetics) method and the classical-PI method are illustrated in 
Figures 5, 6, 7, and 8. 
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Figure 8. The turbine aerodynamic power 


In Figure 5, we can see that the MPPT (PI-Genetics) can quickly converge to MPPT, the power 
factor of MMPT PI-Genetics could finally converge to MMPT at 10ms while the MPPT Pl-classic converges 
towards MPPT at 20ms.It concluded that the Pl-Genetic MPPT could faster track the MPPT than the 
MPPT PI-classic. 

On the other hand, in the context of power factor stabilization, the MPPT PI-classic has suffered a 
significant loss and caused an oscillation around the MPPT, while the MPPT PI-Genetics also offers the 
possibility to eliminate the oscillation around the MPPT. We found that the efficiency of the systemwith 
MPPT PI-Genetics was greater than 99.97%.Based on Figure 5 and Figure 7, it can be seen that PI-Genetics 
is able to reduce the speed and the disturbingpower after the wind turbine has been approached. From the 
obtained simulation results, the MPPT PlI-Genetic has a better performance compared to the MPPT 
Pl-classic. 


6. ELECTRICALCONVERSION CHAIN MODELING 

The chosen Generator for the wind energy transformation is the doubly fed induction generator 
(DFIG) [20, 21].Moreover, a DFIG is controlled by the rotor and is used with a speed variation range 
limitation +50% of the nominal speed. 

The DFIG dynamic model in the reference dq which is related to the rotating field, is presented by 


V V ; i ; 1 
V V i i 
(equation number), where 2° , dr, 2 ds | @ dr | @" are respectively the voltages and the 
currents at the output of the generator in the Park model [22-25]. 


d 
Vas = Rsigs + — Pgs — sb qs 
dt 
d 
Vas = Rigs + — ¢g5 + OP i. 
dt 
d 
Vi = Ryi gy + PAC + (@> — w) ddr 


d 
Vor = Ryigr + Pac (0. = O) 4, 


(15) 
In addition, the flowing equations are expressed by: 
Pas = Lsigs + Mig, 
bg5 = Lsigs + Mig, 
Pay = Lyigs + Migs 
bgr= Lyigr + Migs (16) 


On the other hand, the electromagnetic torque equation is given by: 
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M : 
Cem = > P—— Gastqr — Pqs'ay? 
2. de 


(17) 


The stator flux orientation allows the magnitudes decoupling of the two axes d and q. By orienting 
the stator flux along the axis d, we can write 


sq = 0 Os = Ped 
The electromagnetic torque is then written: 


3 M 
Cem = — P— Gastar? 
2 L, (18) 


In addition, the flows equations become: 


QO=L5 5 t+ Mi 


i = Leig, + Mig, 
q 


= (19) 


In the two-phase reference system, the active and reactive stator powers of an asynchronous 
generator are written: 


Pe Vistds + Vasigs 
O20 iy Vs 
qs ds ds qs (20) 


The adaptation of these equations to the simplifying assumptions gives: 








M . 
P's Ve 7, bn 
Ls 
Vio VM 
O = STS 7 oe 
L Lb 
S S (21) 


For the generator control, the established expressions showing the relationship between the rotor 
currents and voltages that will be applied to it 


Var = Rypigy — 85 (Ly, — 5 Ngr 


S SS (22) 


7. CONVERTERS and ELECTRICITY NETWORK MODEL 
The two converters interposed between the DFIG rotor and the network is the PWM type. 
The control signals will be determined by serving the DC bus voltage [26]. 


V v |v 


aa ae wdl 
ae ee 
Sd wq1 (23) 
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V V 
md | a wd 2 
V 2 |V 
mq wd 2 (24) 


The voltage of the continuous bus u varies according to the power exchanging between the wind 
turbine and the network, it evolves according to the following equation: 


dV 
de _ ly — Tres | 


dt C (25) 





The power that should be circulated by the rotor converter AC / DC / AC in the hyper or hypo 
synchronous case is given by [27]: 


P, 7 ee = § max Ps (26) 


The DC bus voltage must check the following equation: 


V. =2av V2 
dc S (27) 
8. CONTROL of the SENT POWERS to the NETWORK 
The transited active and reactive powers are given as follows: 
‘0 Vardar + Vile 
O = V orhar = Valor (28) 


The currents on the grid side are envisaged by imposing references for the active and 
reactive power, 


V + V 
P ref dr — mes © ef qr — mes 
dr—ref = v2 Fe y2 
dr — mes qr — mes 


V — V 
_ ref qr—-—mes © ef dr — mes 


‘gr —ref — v2 fa y2 
dr — mes qr — mes (29) 
| a 1 i —P 
The reference power: ref ee ee (30) 
P =C V —V 
reg u (sx d — ref ae (31) 


9. SIMULATION RESUTS and DISCUSSION 

All simulations were performed under MATLAB / SIMULINK software environment on a turbine 
system as shown in Figure 1.Two MPPT techniques are implemented in simulation as a controller, to study 
and compare the dynamic response of the turbine system. 
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The reference of the reactive power of 1s respectively set at (SO00; 0; -5000) VAR. The active 
power is fixed at 48% of the mechanical power of the wind turbine. The reference of the DC bus voltage is 
set at ISOV. 

The result of the control of the active power is presented in Figure 8, where it gives perfectly the 
image of the corrected and corrugated currents. The DC bus voltage control result is shown in Figure 10 
where the measured voltage is identical to the reference. The reactive power control result is shown 
in Figure 9. 

Between 0 and 1 s, the system provides the reactive power (Q = SOOOVAR); from | to 3 s, 
the system has the unit power factor (Q = O) and from 3 to 4s the system absorbs the reactive power (Q = - 
5000V AR). 


So 
for) 






































i 
ES a 
ro 
S 
S i) 








Stator power (W) 








°o 




















Classical PI 








(2) 








Classical P| 








Genetical PI 








2.00! 


2.015 








i=) 

















Stator reactive power (w) 

















Genetical P| 
: 0.05 0.4 0.15 0.2 2 


25 3 3.5 4 







































































0 0.5 1.5 





2.5 3 3.5 4 


Time (8) 


Figure 9. The active power of the stator Figure 10. The reactive power of the stator 
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Figure 11. The continuous bus voltage Figure 12. The measured current on the DC bus 


10. CONCLUSION 

This article proposes a new MPPT control for the wind energy generation using the Genetic 
Algorithm approach. To evaluate the effectiveness of this newest control, the MPPT Genetical Pl 
performance was compared to the Classical PI. Finally, the MPPT Genetical PI can more quickly follow the 
MPPT and control the wind system for better performance in terms of power stability. 
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APPENDIX 
Table 1. Wind Energy Conversion System Parameters Used 
Parameter Nominal value 
Ray R 35 1 
total inertia J 0.3125 kg.m2 
Coefficient viscous f, 0.00681 
Multiplier G 5.4 
Power P,, 7.5 KW 
Effective stator voltage V,, 220V 
Stator frequency f 50Hz 
stator resistor R, 0.455 0 
rotor resistor R, 0.62 O 
magnetizing inductance 0.078 H 
stator inductors 0.048 H 
rotor inductors 0.081 H 
population size 50 
Number of variables 2 
probability of mutation 5% 
Number of generation 50 
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